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Abstract 


Work  during  this  period  was  concerned  with  two  aspects:  1)  development  of  a  new 
method  for  the  derivation  of  the  state  equations  of  motion  for  control  of  flexible 
spacecraft  and  2)  development  of  a  method  for  the  control  of  spacecraft  in  the  form  of 
articulated  flexible  multi-bodies. 

In  deriving  the  equations  of  motion  for  the  control  of  flexible  spacecraft,  it  is  common 
practice  to  express  the  rigid-body  rotational  motions  in  terms  of  nonorthogonal  inertial 
components,  such  as  Euler’s  angles.  On  the  other  hand,  in  controlling  the  same  motions 
of  the  spacecraft,  the  control  law  is  generally  in  terms  of  angular  motions  components 
along  the  orthogonal  body  axes.  The  elastic  motions  are  ordinarily  expressed  in  terms 
of  components  along  body  axes.  In  implementing  the  control  laws,  it  is  desirable  that 
both  the  equations  of  motion  and  the  control  laws  be  in  terms  of  the  same  type  of  var¬ 
iables,  and  in  particular  in  terms  of  the  body  axes  variables.  Reference  1  presents  the 
development  of  state  equations  of  motion  for  flexible  bodies  in  terms  of  quasi¬ 
coordinates,  which  represent  this  very  type  of  variables.  Indeed,  the  angular  velocity 
components  about  the  orthogonal  body  axes  are  time  derivatives  of  quasi-coordinates. 
The  state  equations  of  motion  for  flexible  bodies  in  terms  of  quasi-coordinates  render 
the  control  implementation  task  considerably  easier. 

In  many  space  applications,  it  becomes  necessary  to  reorient  the  line  of  sight.  If  the  line 
of  sight  is  fixed  in  the  spacecraft,  such  as  in  the  space  telescope,  this  implies  reorien¬ 
tation  of  the  whole  spacecraft.  In  other  space  applications,  the  spacecraft  consists  of  a 
rigid  platform  and  a  number  of  flexible  appendages,  such  as  antennas.  If  the  mission 
requires  the  reorientation  of  the  line  of  sight  of  these  antennas  relative  to  the  inertial 
space,  then  the  preferred  strategy  is  to  stabilize  the  platform  relative  to  the  inertial  space 
and  to  reorient  the  flexible  antennas  relative  to  the  platform.  Reference  2  provides  a 
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derivation  of  the  state  equations  of  motion  for  the  task  described  above,  based  on  the 
equations  of  motion  in  terms  of  quasi-coordinates  derived  in  Ref.  I.  Because  the 
motions  defining  the  maneuvering  of  the  antennas  arc  known  a  priori,  the  state 
equations  contain  time-dependent  coefficients  and  persistent  disturbances. 

The  problem  of  maneuvering  an  articulated  flexible  spacecraft  and  controlling  its  vi¬ 
bration  at  the  same  time  is  a  very  complex  task.  Reference  3  shows  an  approach  to  the 
problem,  in  which  the  maneuvering  of  the  appendages  is  carried  out  open-loop  using  a 
bang-bang  control  law.  On  the  other  hand,  the  vibration  suppression  is  carried  out 
closed-loop,  which  amounts  to  controlling  a  time-varying  system  subjected  to  persistent 
disturbances. 

In  addition  to  the  above  research,  work  was  carried  out  on  several  other  papers  (Refs. 
4  -  7)  describing  research  under  the  preceding  AFOSR  grant.  In  particular,  we  single 
out  Ref.  5,  which  was  presented  during  this  grant  period.  The  paper  is  concerned  with 
control  of  the  perturbations  experienced  by  a  flexible  spacecraft  during  a  minimum-time 
maneuver.  The  spacecraft  is  modeled  as  a  flexible  appendage  attached  to  a  rigid  hub. 
The  perturbed  model  can  be  divided  into  a  rigid-body  part  and  an  clastic  part.  The 
model  is  described  by  a  linear,  time-varying  set  of  ordinary  differential  equations  sub¬ 
jected  to  piecewise-constant  disturbances  caused  by  inertial  forces  resulting  from  the 
minimum-time  maneuver.  The  control  consists  of  a  reduced-order  compensator  de¬ 
signed  such  that  the  perturbed  model  is  finite-time  stable. 
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State  Equations  of  Motion  for  Flexible  Bodies 
in  Terms  of  Quasi-Coordinates 
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Summary 

This  paper  is  concerned  with  the  general  motion  of  a  flexible 
body  in  space.  Using  the  extended  Hamilton's  principle  for  dis¬ 
tributed  systems,  standard  Lagrange's  equations  for  hybrid  sys¬ 
tems  are  first  derived.  Then,  the  equations  for  the  rigid-body 
motions  are  transformed  into  a  symbolic  vector  form  of  Lagrange’s 
equations  in  terms  of  general  quasi-coordinates.  The  hybrid 
Lagrange's  equations  of  motion  in  terms  of  general  quasi-coordi¬ 
nates  are  subsequently  expressed  in  terms  of  quasi-coordinates 
representing  rigid-body  motions.  Finally,  the  second-order 
Lagrange's  equations  for  hybrid  systems  are  transformed  into  a 
set  of  state  equations  suitable  for  control.  An  illustrative 
example  is  presented. 


Introduction 

The  derivation  of  the  equations  of  motion  has  preoccupied  dynam- 
icists  for  many  years,  as  can  be  concluded  from  the  texts  by 
Whittaker  [1],  Pars  [2]  and  Meirovitch  [3].  References  1-3  con¬ 
sider  the  motion  of  systems  of  particles  and  rigid  bodies,  and 
the  equations  of  motion  are  presented  in  a  large  variety  of 
forms.  In  this  paper,  we  concentrate  on  a  certain  formulation, 
namely,  Lagrange's  equations.  For  an  n-degree-of-f reedom  system, 
Lagrange's  equations  consist  of  n  second-order  ordinary  differen¬ 
tial  equations  for  the  system  displacements. 

In  the  control  of  dynamical  systems,  it  is  often  convenient  to 
work  with  first-order  rather  than  second-order  differential  equa- 
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tions.  Introducing  the  velocities  a3  auxiliary  variables,  it  is 
possible  to  transform  the  n  second-order  equations  into  2n  first- 
order  state  equations.  The  state  equations  are  widely  used  in 
modern  control  theory  [4]. 

With  the  advent  of  man-made  satellites,  there  has  been  a  renewed 
interest  in  the  derivation  of  the  equations  of  motion.  The 
motion  of  rigid  spacecraft  can  be  defined  in  terms  of  transla¬ 
tions  and  rotations  of  a  reference  set  of  axes  embedded  in  the 
body  and  known  as  body  axes.  The  equations  of  motion  for  such 
systems  can  be  obtained  with  ease  by  means  of  Lagrange's  equa¬ 
tions.  It  is  common  practice  to  define  the  orientation  of  the 
body  relative  to  an  inertial  space  in  terms  of  a  set  of  rotations 
about  nonorthogonal  axes  [3].  However,  the  kinetic  energy  has  a 
simpler  form  when  expressed  in  terms  of  angular  velocity  compo¬ 
nents  about  the  orthogonal  body  axes  than  in  terms  of  angular 
velocities  about  nonorthogonal  axes.  Moreover,  for  feedback 
control,  it  is  more  convenient  to  work  with  angular  velocity 
components  about  the  body  axes,  as  sensors  measure  angular 
motions  and  actuators  apply  torques  in  terms  of  components  about 
the  body  axes.  In  such  cases,  it  is  often  advantageous  to  work 
not  with  standard  Lagrange's  equations  but  with  Lagrange's  equa¬ 
tions  in  terms  of  quasi-coordinates  {1,3].  If  the  body  contains 
discrete  parts,  such  as  lumped  masses  connected  to  a  main  rigid 
body  by  massless  springs,  it  is  convenient  to  work  with  a  set  of 
axes  embedded  in  the  undeformed  body.  The  equations  of  motion 
consist  entirely  of  ordinary  differential  equations  and  can  be 
obtained  by  a  variety  of  approaches,  including  the  standard 
Lagrange's  equations  and  Lagrange's  equations  in  terms  of  quasi¬ 
coordinates  [5]*. 

In  the  more  general  case,  the  body  can  be  regarded  as  being 
either  entirely  flexible  with  distributed  mass  and  stiffness 
properties  or  as  consisting  of  a  main  rigid  body  with  distributed 
elastic  appendages.  Unlike  the  previous  case,  the  equations  of 
motion  are  hybrid,  in  the  sense  that  the  equations  for  the  rigid- 


*  Note  that  Ref.  5  refers  to  Lagrange's  equations  in  terms  of 
quasi-coordinates  as  Boitzmann-Hamel  equations. 
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body  motions  are  ordinary  differential  equations  and  those  for 
the  elastic  motions  are  partial  differential  equations.  Hybrid 
equations  were  obtained  for  the  first  time  in  Ref.  6.  Moreover, 
the  formulation  of  Ref.  6  was  obtained  by  using  Lagrange's  equa¬ 
tions  in  terms  of  quasi-coordinates,  but  some  generality  was  lost 
in  that  the  body  considered  was  assumed  to  be  symmetric  and  to 
undergo  antisymmetric  elastic  motion.  As  a  result,  the  rigid- 
body  translations  were  zero. 

This  paper  is  concerned  with  the  general  motion  of  a  flexible 
body  in  space.  Using  the  extended  Hamilton's  principle  for  dis¬ 
tributed  systems  (7),  standard  Lagrange's  equations  for  hybrid 
systems  are  first  derived.  Then,  using  the  approach  of  Ref.  3, 
the  equations  for  the  rigid-body  motions  are  transformed  into  a 
symbolic  vector  form  of  Lagrange's  equations  in  terms  of  general 
quasi-coordinates.  The  hybrid  Lagrange's  equations  of  motion  in 
terms  of  general  quasi-coordinates  are  subsequently  expressed  in 
terms  of  quasi-coordinates  representing  rigid-body  motions.  This- 
is  a  very  important  step,  as  the  latter  form  permits  the  derivat¬ 
ion  of  the  hybrid  equations  of  motion  with  relative  ease,  thus 
eliminating  a  great  deal  of  tedious  work.  These  hybrid  equations 
represent  an  extension  to  flexible  bodies  of  Lagrange's  differen¬ 
tial  equations  in  terms  of  quasi-coordinates  derived  in  Ref.  3 
for  rigid  bodies.  The  second-order  equations  are  then  used  to 
derive  the  hybrid  state  equations. 

As  an  illustration,  the  hybrid  equations  of  motion  of  a  space¬ 
craft  consisting  of  a  rigid  hub  with  a  flexible  appendage  simu¬ 
lating  an  antenna  are  derived. 

Standard  Lagrange’s  Equations  for  Hybrid  Systems 

Let  us  consider  a  flexible  body  and  assume  that  the  Lagrangian  L 
=  T  -  V,  in  which  T  is  the  kinetic  energy  and  V  is  the  potential 
energy,  can  be  written  in  the  general  form  L  -  L  (q  .  ,q  .  ,u  .,u  ;,u‘.  ,u‘! , . . 
..Uj  ).  where  qj  =  q  ^  ( t )  (i  =  1,2,.  ...m)  are  generalized  coordinates 
describing  rigid-body  motions  of  the  body  and  Uj(P,t)  (j  = 

1,2,. ..,n)  are  generalized  coordinates  describing  elastic  motions 
relative  to  the  rigid-body  motions  of  a  typical  point  in  the  body 
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identified  by  the  spatial  position  P.  Dots  designate  derivatives 
with  respect  to  time  and  primes  derivatives  with  respect  to  the 
spatial  position.  for  convenience,  we  express  the  Lagrangian  in 
terms  of  the  Lagrangian  density  L  in  the  form  L  =  j\,  L  dO,  where  0 
is  the  domain  of  extension  of  the  body. 


We  propose  to  derive  Lagrange's  equations  by  means  of  the 
extended  Hamilton's  principle  [7],  which  can  be  stated  as 


fD  («L  +  «W)dDdt  =  0.  5q i  =  dUj  =  0  at 


tl*  t2 


(1) 


where  «W  is 
related  to 
be  written 


the  nonconservative  virtual  work 
the  virtual  work  by  6W  =  Jg  5W  dO. 
in  the  form 


density,  which  is 
The  virtual  work  can 


6U 


L  Vq1  +  l  .'o  Uj{uj  d0 

i=l  j=l  J 


(2) 


where  are  nonconservative  generalized  forces  associated  with 
the  rigid  body  motions  and  are  nonconservative  generalized 
force  densities  associated  with  the  elastic  motions;  aq^  and  5Uj 
are  associated  virtual  displacements.  Following  the  usual  steps 
[7),  we  obtain  Lagrange's  equations  of  motion,  which  can  be 
expressed  in  the  symbolic  vector  form 


a_  f al 

V  ' 


Q. 


i_  fit)  _ 

at  v  • 1 

3U 


il 

3U 


£u  =  U 


( 3a,b) 


where  g  and  Q  are  m-vectors,  u  and  0  are  n-vectors  and  L  is  an 
n  «  n  operator  matrix.  Because  of  the  mixed  nature  of  the  differ¬ 
ential  equations,  we  refer  to  the  set  (3)  as  hybrid.  The  elastic 
displacements  are  subject  to  given  boundary  conditions. 


Equations  in  Terms  of  Quasi-Coordinates  for  the  Rigid-Body 
Motions 

Quite  often  it  is  convenient  to  express  the  Lagrangian  not  in 
terms  of  the  velocities  q^  but  in  terms  of  linear  combinations 
(t=l,2, ...  ,m)  of  qf.  The  difference  between  q^  and  Wj  is  that  the 
former  represent  time  derivatives  dq^/dt,  which  can  be  integrated 
with  respect  to  time  to  obtain  the  displacements  q^,  whereas 

cannot  be  integrated  to  obtain  displacements.  It  is  customary 
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to  refer  to  as  derivatives  of  quasi-coordinates  [3],  The 
relation  between  Wg  and  q^  can  be  expressed  in  the  compact  matrix 
form  w  «  A^q,  where  the  notation  is  obvious.  Similarly,  we 
express  the  velocities  q^  in  terms  of  the  variables  w  as  5  =  6w, 
from  which  it  follows  that  the  m  «  m  matrices  A  and  B  are  related 
by  A  8  =  B^A  =>  I,  where  I  is  the  identity  matrix  of  order  m. 


Our  object  is  to  derive  Lagrange's  equations  in  terms  of  w 
instead  of  q^.  Using  the  relations  indicated  above,  it  can  be 
shown  [3]  that  Eqs.  (3a)  can  be  replaced  by 


d_  /3L,  „TP  3L  „T  3L  _  M 
dt  l  wJ  +  8  E  aS  -  8  3q  '  - 


where 


|/bt 


l-V  U-|.  N  -  BTQ 


(4) 


(5a, b) 


and  we  note  that  the  first  matrix  in  £  is  obtained  by  first 
carrying  out  a  triple  matrix  product  for  every  one  of  the  m 
entries  in  A  and  then  arranging  the  resulting  scalars  in  a  square 
matrix.  On  the  other  hand,  the  second  matrix  in  E  is  obtained  by 
first  generating  a  row  matrix  for  every  generalized  coordinate  q^ 
(k  *  1,2,... ,m)  and  then  arranging  the  row  matrices  in  a  square 
matrix.  Equation  (4)  represents  a  symbolic  vector  form  of  the 
Langrange  equations  for  quasi-coordinates.  The  complete  formula¬ 
tion  is  obtained  by  adjoining  to  Eq.  (4),  the  equations  for  the 
elastic  motion,  Eq.  (3b),  as  well  as  the  associated  boundary 
conditions . 


General  Equations  in  Terms  of  Quasi-Coordinates  for  a 
Translating  and  Rotating  Flexible  Body. 

Let  us  consider  the  body  depicted  in  Fig.  1.  The  motion  of  the 
body  can  be  described  by  attaching  a  set  of  body  axes  xyz  to  the 
body  in  undeformed  state.  The  origin  of  the  body  axes  coincides 
with  an  arbitrary  point  0.  Then,  the  motion  can  be  defined  in 
terms  of  the  translation  of  point  0,  and  the  rotation  of  the  body 
axes  xyz  relative  to  the  inertial  axes  XYZ.  The  position  of  0 
relative  to  XYZ  is  given  by  the  radius  vector  R  =  R(RX,  Ry,  R2). 

The  rotation  can  be  defined  in  terms  of  a  set  of  angles  8^,  3., 

and  8j  (Fig.  2).  Hence,  the  generalized  coordinates  are 

^1  *  By*  Rj  *  By*  R3  *  Rj*  44  1  -  92*  ^6  =  e2'  In  Edition,  there 
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ace  the  elastic  displacement  components  u  (P,t),  u  (P,t),  u7(P,t). 

a  y  4 

The  displacements  Ry,  Ry,  R7  are  measured  relative  to  the  inertial 
axes  XYZ.  On  the  other  hand,  the  displacements  u  ,  u.  ,  u7  are  mea- 
sured  relative  to  the  body  axes  *yz.  Moreover,  the  components 
R)(.  Ry.  Rz  of  the  velocity  vector  S  are  also  measured  relative  to 
XYZ.  On  the  other  hand,  the  angular  velocity  vector  j  has  compo¬ 
nents  jy,  measured  'elative  to  the  body  axes  <yz.  It  will 
prove  convenient  to  express  all  motions  in  terms  of  components 
along  the  body  axes.  To  this  end,  if  we  denote  the  velocity  of 
point  0  in  terms  of  components  along  the  body  axes  by  V,  then  it 
can  be  shown  that  V  =  CR,  where  C  -  is  a  rotation  matrix. 

Moreover,  the  angular  velocity  vector  u  can  be  expressed  in  terms 
of  the  angular  velocities  8^,  8^  and  in  the  form  *  =  Oa,  where 
0  *  0(8. ,9,)  is  a  transformation  matrix.  We  note  that  the  angular 

velocity  components  u  ,  and  u,  cannot  be  integrated  with  respect 

x  y  z 

to  time  to  yield  angular  displacements  1  ,  a  and  a  about  axes  x,  y 

x  y  z 

and  z,  respectively.  Hence,  ^  can  be  regarded  as  time 

x  y  4 

derivatives  of  quasi-coordinates  and  treated  by  the  procedure 
presented  in  the  preceding  section.  Although  it  is  not  very 
common  to  regard  the  velocity  components  Vx,  Vy  and  Vz  as  time 
derivatives  of  quasi-coordinates,  they  can  still  be  treated  as 
such.  In  view  of  this,  if  we  introduce  the  generalized  velocity 
vector  q  =  (Ry  Ry  R^e^  a^  8j|\  as  well  as  the  "quasi-velocity" 
vector  w  =  [Vx  V2  wx  u2)\  we  conclude  that  the  coefficient 
matrices  are  defined  by 


C_  j_0 
0  !  0 


BT  *  A"1 


lo 


(oY1 


(6a, b) 


where  we  recognized  that  C'*  =  CT,  because  rotation  matrices 
are  orthonormal.  It  can  be  shown,  after  lengthy  algebraic  manip¬ 
ulations,  that 


8TE  = 


v  i : 


(7) 


where  w  and  V  are  skew-symmetric  matrices  corresponding  to 
j  and  V  (3),  respectively. 


Using  Eqs.  (3b)  and  (4)  in  conjunction  with  the  above  relations, 
we  obtain  the  hybrid  Lagrange's  equations  in  terms  af  quasi-coor 
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dinates 


1-  fii) 

dt  ‘ 3  V 


—  -  C  —  =  F 
;V  L  ifi  - 


-  (-] 

dt  Mu' 


,  4  +  :  ii  .  (0T)-i  it 

3V  -,U  v  '  33 


-  !-) 
3 1  MvJ 


ll 

3U 


+  £  u  =  u 


(  8a) 

(8b) 

(  8c) 


where  F  and  H  are  external  nonconservative  force  and  torque, 
respectively,  in  terms  of  components  aoout  the  body  axes, 

3L/38  -  (3L/38j  3 L/ a e ^  3 L/ 3 e ^ | ^  and  v  =  u.  Note  that  a  does  not 

really  represent  a  vector  and  must  be  interpreted  as  a  mere  sym¬ 
bolic  notation.  We  recall  that  the  components  of  u  are  still 
subject  to  given  boundary  conditions. 


It  should  be  pointed  out  that,  in  deriving  Eqs.  (8),  no  explicit 
use  was  made  of  the  angles  8,,  3^  and  a3>  so  that  Eqs.  (8)  are 
valid  for  any  set  of  angles  describing  the  rotation  of  the  body 
axes,  such  as  Euler's  angles,  and  they  are  not  restricted  to  the 
angles  used  here.  Moreover,  point  0  is  an  arbitrary  point,  not 
necessarily  the  mass  center  of  the  undeformed  body,  and  axes  xyz 
are  not  necessarily  principal  axes  of  the  undeformed  body.  Clear¬ 
ly,  if  xyz  are  chosen  as  the  principal  axes  with  the  origin  at  the 
mass  center,  then  the  equations  of  motion  can  be  simplified. 


State  Equations  in  Terms  of  Quasi-Coordinates 

Equations  (8),  and  in  particular  Eqs.  (8a)  and  (8b),  can  be 
expressed  in  more  detailed  form.  To  this  end,  we  write  the 
velocity  vector  of  a  typical  point  P  in  the  body  in  terms  of  com¬ 
ponents  along  the  body  axes  as  follows: 

Vp  =  V  +  Ux(r  +  u)+v  =  Vt(r  +  u)^u  v  (9) 

where  r  is  the  nominal  position  of  P  relative  to  0.  Moreover, 
r  and  u  represent  skew-symmetric  matrices  associated  with  the 
vectors  r  and  u,  respectively.  Then,  denoting  by  o  the  mass  den¬ 
sity,  the  kinetic  energy  can  be  shown  to  have  the  expression 

T  =  |  J"g  °,p^p  dO  *  |  *  vV.  +  VTf0  ov  dD 


(10) 


+ 


o(r  *  ii)v  dD 


1  T 

2  ? 


Ju  +  ^  Jq  ov^v  dO 


where  S  *  J"q  p(r  +  u ) dO ,  J  =  ,g  p(r  +  u)(r  +  u)^dD,  in  which  S  is 

recognized  as  a  skew-symmetric  matrix  of  first  moments  and  J  as  a 
symmetric  matrix  of  mass  moments  of  inertia,  both  corresponding 
to  the  deformed  body.  Moreover,  we  assume  that  the  potential 
energy  has  the  functional  form  V  *  V(R,  e,  u,  u1 . u^). 


Inserting  Eq.  (10)  into  Eqs.  (8)  and  rearranging,  we  obtain  the 
explicit  Lagrange's  equations  in  terms  of  hybrid  coordinates 

mV  +  STw  +  f0flv  dO  =  (2Sy  ♦  mV  +  iS)u  -  C  +  F  (11a) 

SV  +  Jw  +  Jjj  p(r  +  u)v  dD  *  [2  o(r  +  u)v  dD  +  SV  -  uJju 

-  (O1)"1  If  ♦  ?  (lib) 

-  -  T*  *  "I  -2  -T 

oV  +  o(r  +  u)  u  +  pv  =  -  3V  ui  -  au  (r  +  u)  -  2dv  u  -  Lu  +  U  (11c) 


where  S  *  f  pv  dO.  The  state  equations  are  completed  by 
D 

adjoining  the  kinematical  relations 


CTV,  8 


0  , 


U  =  V 


( lid, e, f  ) 


Illustrative  Example 

As  an  illustration,  we  consider  a  spacecraft  consisting  of  a 
rigid  hub  and  a  flexible  appendage,  as  shown  in  Fig.  3.  From 
the  figure,  we  can  write 


r  -  xi, 
so  that 


u  =>  uyj  +  uzk,  v  =  vyJ  ♦  vzk 


(12) 


J"puzdx 

-fouydx 


-/ouzdx 

[puzdx 


mx 


fpuydx' 

-mx 

0 


(13a) 


where  o  is  the  mass  density  of  the  appendage,  m  is  the  total  mass  and  x  is  the 
position  of  the  mass  center  of  the  appendage.  Moreover, 
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"jxx+rp(u^u^dx 

J  =  -foxu  dx 
y 

-Joxu^dx 


-foxu^dx  -fpxi^dx 

O 

Jyy^oUzdx  -r»Vzdx 

-;Du  u2dx  JZz^‘,uydx. 


where  ,  J  and  Jzz  are  the  mass  moments  of  inertia  of  the 
spacecraft  regarded  as  rigid. 

Using  Eqs.  (12)  and  (13),  the  state  equations,  Eqs.  (11),  can  be 
written  in  the  explicit  forms 

Rx  «  (ce2c93  +  S91se2se3)Vx  -  (cSjCa^  -  sSjSSjCe^V  +  ce^sa^  (14a) 


Ry  *  c9iS03vx  *  cejCSjV  -  sa^ 


(14b) 


Rz  *  -(S92C93  -  59^92592)7^  +  (S82S93  +  +  08^9^  (14c) 


.  .  59  3  C8  3 

9 1  =  C9-J  tu  -  S90  u,  i  6a  s  • —  w  +■  ■  ~ —  (j 

1  3  x  3  y*  2  C9|  x  C9^  y 

.  s0is03  se.ce^ 

3  C8,  x  C9,  y  z*  y  y’  z  z 


( 14d,e) 


( 14f ,g,h) 


0  0 

mV^  +  u>yJpU2  dx  -  dx  =  mVyWz  *  m^zuy  +  n,ix(“y  *  uz)  ”  ux“yJou 


-  <V2foUz  dx  +  2uzJ0vy  dx  -  2wxfpvz  dx  -  (c82C9j  +  59^62583) 


-  C92S93  %  ♦  <se2ce3  -  S9JC92S9;,)  fgj  ♦  Fx 


( 1 4  i ) 


—  _  o  0 

mVy  -  <*>xfpuz  dx  +  m3xuz  *  m^zux  ■  *>VX«Z  -  iti^xu^u^  +  (wx  +  “z)Jpu^  dx 


-  <y.z/p U2  dx  *  2uxJpVz  dX  ♦  (C92S83  -  S9lS92C93)  |jj- 


-  (s92S93  +  S8lC92C93)  w  +  F 


(14j) 


—  1— 1  0  0 

mVz  +  “x^°uy  dx  *  mlXuy  =  mVxuy  ■  mVywx  *  mlXux“z  +  (“x  +  uyIoUz  dx 


r  r  g  v  iu  ay 

-  Vzf0Uy  dx  *  2uxIoVy  dx  -  C9lS92  ^  -  C9lC92  -fq  *  Fz 


(14k) 
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(r0uz  dx)Vy  +  (Jouy  dx)Vz  +■  (Jxx  *  fo(uy  t-  Uz)dxlux  -  (|0xuy  dx)J,y 

-  (Joxuz  dx)uz  *  fo(uyvz  -  uzvy)dx  =  -  (Vy.foUy  dx  ♦  vzl^uz  dx) 

+  V  111  foil  dx  +  V  -j  oil  dx  +  j  J  oXll  dx  -  u  a  "jXll  dx 
x  yJ  y  x  z-  z  *  y  *  z  a  Z"  y 

+  (u.y  -  *z)>yuz  dx  +  JyJzIJyy  -  Jzz  -  ,'o(uy  -  u2)dx| 

se,  3 V  sa.se,  u 

-  2u  fo(u  v  t  u  v  )dx  -  C3,  — - - — - - — -  - —  <-  M  (141) 

XJ  '  y  y  2  2J  3  33  J  C8j  33^  C91  "3  x 


_  •  n  •  • 

(JoUzdx)Vx  -  m1xVz  -  (j'oxuz  dx)cjx  *  (Jyy  +  /ouz  dx)Uy  -  (/o uyuz  dx)u,? 


J2Vx  *  xv2)dx  5  VV-1*  -  (VZJ*°‘JZ  dx  +  m1xVx)u  +  V  oz|ouz  dx 


-  “7"«(d*v  -  d77  +  loU£  dx)  -  (<*>*  -  u7 ) j*o xii  dx  +  u  u Joxu  dx 


Z  Xx  XX  zz 


x  z/J  z  y  zJ  y 


-  “xwy^0UyU2  dx  +  2lJx-roXVy  dx  '  2ayi"(3U2v2  dx 

.  3W  c®3  au  S6iC3,  y 

+  2u  foU  V  dx  +  S3,  —  -  - -  - - - - -  —  +  M 

zJ  z  y  3  ae1  ceL  se2  ca1  aa^  y 


(14m) 


-  (IoUy  dx)Vx  +  m1xVy  -  (Joxuz  dx)ux  -  (fouyuz  dx)Uy  +  (Jzz  +  f0uy  dx)U)< 
+  J‘o(xvy  -  uyvx)dx  =  m1xVzux  +  VzuyJouy  dx  -  (VyJouy  dx  -  V^x)^ 

p  p  p 

*  “x“y^Jxx  “  Jyy  +  ^olJy  dx^  +  (“x  *  “y)I°XUy  dx  *  “yu2-'3Xuz  dx 


+  v.z/oUyuz  dx  +  2u.x/oxvz  dx  +  2Uy/0uyvz  dx  -  2«z/0uyvy  dx  - 


(14n) 


•  2  2 
oVy  -  oUzwx  +  °Xu,z  *  °vy  *  °^z<4,x  ”  °^xwz  "  pXwxwy  +  °(wx  *  wz^uy 


-  o«yUzuz  +  2ovz«x  -  4yuy+  uy 


oVZ  +  oUyux  *  oX“y  +  °v2  *  *  aVylJx  >  oVx“y  *  oX“xli,z  '  ouy“zuy 

♦  o(«x  +  wy)uz  -  2ovyux  -  £zuz  +  Uz  (14p) 

where  is  the  mass  o£  the  appendage,  se  ^  =  sin  e  ^ ,  ce.  =  cos  a  ^  ( i  => 

1,2,3)  and 
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of  elasticity  and 

and  I2  are  area 

-  '-2  -  77 

3  X  ix 


moments  of  inertia.  The  operators  L y  and  £~  z  include  the 
effects  of  bending  and  of  the  axial  force  on  the  appendage  [7' 


Summary  and  Conclusions 

In  deriving  the  equations  of  motion  for  flexible  bodies  by  the 
Lagrangian  approach,  it  is  common  practice  to  express  the  rota¬ 
tional  motion  in  terms  of  angular  velocities  about  nonor t hogona 1 
axes,  which  tends  to  complicate  the  equations.  Moreover,  this 
creates  difficulties  in  feedback  control,  in  which  the  torque 
actuators  apply  moments  about  body  axes  and  the  output  of  sensors 
measuring  angular  motion  is  also  expressed  in  terms  of  components 
about  the  body  axes.  The  same  can  be  said  about  force  actuators 
and  translational  motion  sensors.  It  turns  out  that  the  equat¬ 
ions  of  motion  are  appreciably  simpler  when  the  rigid-body  trans¬ 
lations  and  rotations  are  expressed  in  terms  of  components  about 
the  body  axes.  Such  equations  can  be  obtained  by  introducing  the 
concept  of  quasi-coordinates.  The  concept  of  quasi-coordinates 
was  used  earlier  by  this  author  to  derive  equations  of  motion  of 
rotating  bodies  with  flexible  appendages,  but  never  in  the 
general  context  considered  here.  Indeed,  in  this  paper, 
Lagrange's  equations  in  terms  of  quasi-coordinates  are  derived 
for  a  distributed  flexible  body  undergoing  arbitrary  rigid-body 
translations  and  rotations,  in  addition  to  elastic  deforma¬ 
tions.  The  second-order  differential  equations  in  time  for  the 
hybrid  system  are  then  transformed  into  a  set  of  hybrid  state 
equations  suitable  for  control  design.  The  approach  is  demon¬ 
strated  by  deriving  the  hybrid  state  equations  of  motion  for  a 
spacecraft  consisting  of  a  rigid  body  with  a  flexible  appendage 
in  the  form  of  a  beam. 
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Figure  3.  A  Rigid  Spacecraft  with  a  flexible  Appendage 
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This  paper  is  concerned  with  the  derivation  of  the  state  equations  of  motion  for 
a  spacecraft  consisting  of  a  main  rigid  platform  and  a  given  number  of  flexible 
appendages  changing  the  orientation  relative  to  the  main  body.  The  equations 
are  derived  by  means  of  I.agrange's  equations  in  terms  of  quasi-coordinates. 
Assuming  that  the  appendages  represent  distributed-parameter  members,  the 
state  equations  of  motion  are  hybrid.  Moreover,  they  are  nonlinear.  Following 
spatial  discretization  and  truncation,  the  hybrid  equations  reduce  to  a  system 
of  nonlinear  discretized  state  equations,  which  are  more  practical  for  numerical 
calculations  and  control  design.  To  illustrate  the  effect  of  nonlinearity  on  the 
dynamic  response  during  reorientation,  a  numerical  example  involving 
spacecraft  with  a  membrane-like  antenna  is  presented. 


I.  INTRODUCTION 

In  many  space  applications,  it  becomes  necessary  to  reorient 
a  certain  line  of  sight  in  a  spacecraft.  Examples  of  this  are  the 
reorientation  of  a  space  telescope  or  of  an  antenna  in  a 
spacecraft.  In  some  cases,  such  as  in  the  space  telescope,  the 
line  of  sight  can  be  regarded  as  being  fixed  relative  to  the 
undeformed  structure,  in  which  cases  reorientation  of  the  line 
of  sight  implies  maneuvering  of  the  whole  spacecraft  [Turner 
and  Junkins  (1980),  Dreakwell  (1981),  Raruh  and  Silverberg 
(1984),  Meirovitch  and  Quinn  (1987),  Meirovitch  and  Sharony 
(1987)  and  Quinn  and  Meirovitch  (19S8)|.  However,  many 
spacecraft  can  be  represented  by  mathematical  models 
consisting  of  a  rigid  platform  with  one  or  more  flexible 
appendages,  such  as  flexible  antennas,  so  that  the  mission 
involves  the  maneuvering  of  a  hybrid  (lumped  and  distributed 
flexible)  system.  Quite  often,  the  line  of  sight  coincides  with 
an  axis  fixed  in  a  small  component  of  the  spacecraft,  such  as 
an  antenna,  in  which  case  it  may  be  more  advisable  to  retarget 
only  the  antenna  and  not  the  entire  spacecraft.  This  is 
particularly  true  when  the  inertia  of  the  antenna  is  much 
smaller  than  the  inertia  of  the  spacecraft.  The  argument 
becomes  even  stronger  when  several  antennas  must  be 
retargeted  independently  in  space.  In  such  cases,  it  appears 
more  sensible  to  conceive  of  a  spacecraft  consisting  of  a 
platform  stabilized  in  an  inertial  space  with  several 
appendages,  rigid  or  flexible,  hinged  to  the  platform  and 
capable  of  pivoting  about  two  orthogonal  axes  relative  to  the 
platform  [Meirovitch  and  Kwak  (1988a.  1988b)  and 
Meirovitch  and  France  (1989)|.  In  this  case,  reorientation 
relative  to  the  stabilized  platform  is  equivalent  to  retargeting  in 
an  inertial  space.  Note  that  such  a  maneuvering  spacecraft  is 
characterized  by  the  fact  that  its  configuration  varies  with  time. 
This  paper  is  concerned  with  the  mission  of  independent 
retargeting  of  the  line  of  sight  of  each  antenna  relative  to  the 
inertiai  space. 

Otic  problem  encountered  in  the  dynamics  of  a  complex 
system  is  how  to  derive  the  equations  of  motion  efficiently.  In 
general,  the  equations  of  motion  for  a  flexible  spacecraft  have 


very  complicated  expressions  [Grote,  McMunn  and  Gluck 
(1970)  and  Likins  (1972)).  so  that  new  methods  for  deriving 
equations  of  motion  have  been  proposed  [Ho  (1977)  and  Kane 
and  Levinson  (1980)|.  Kane  and  Levinson  compared  seven 
different  methods.  Lagrange's  equations  of  motion  in  terms 
of  quasi -coordinates  for  a  hybrid  system  were  derived  first  by 
Meirovitch  (196b)  and  then  by  Williams  (1976)  and  Brown 
(1981).  Recently.  Meirovitch  (1988)  and  Meirovitch  and 
Kwak  (1989)  showed  that  Lagrange's  equations  of  motion  in 
terms  of  quasi  coordinates  are  quite  useful  for  deriving  the 
equations  of  motion  for  the  maneuvering  and  control  of 
flexible  spacecraft.  Recause  the  derived  equations  of  motion 
are  based  on  body-fixed  coordinates,  control  design  based  on 
such  equations  is  very  convenient. 

The  mathematical  formulation  for  a  spacecraft  including  a 
rigid  platform  and  several  flexible  antennas  consists  of  a  hybrid 
set  of  equations  of  motion,  in  the  sense  that  there  are  six 
ordinary  differential  equations  for  the  rigid -body  translations 
and  rotations  of  the  platform  and  partial  differential  equations 
for  the  elastic  motion  of  each  antenna.  The  equations  of 
motion  are  not  only  hybrid,  but  the  maneuvering  of  the 
antennas  relative  to  the  platform  according  to  some  prescribed 
function  of  time  introduces  time-dependent  coefficients  into  the 
equations.  Moreover,  the  equations  contain  terms  reflecting 
persistent  disturbances  caused  by  inertial  forces.  Because  both 
numerical  simulation  and  control  design  of  systems  governed 
by  sets  of  hybrid  differential  equations  are  not  feasible,  it  is 
necessary  to  discretize  the  partial  differential  equations  in 
space,  which  can  be  carried  out  by  the  classical  Raylcigh-Rilz 
method  or  the  finite  element  method  [Meirovitch  (1980)). 


II.  GENERAL  LACRANGE'S  EQUATIONS  IN 
TERMS  OF  QUASI-COORDINATES 

L  et  us  consider  a  system  consisting  of  a  main  rigid  body, 
acting  as  a  platform,  and  a  certain  number  of  flexible 
appendages  binged  to  the  main  rigid  body.  The  interest  lies  in 
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reorienting  the  flexible  appendages  independently  so  as  to  point 
in  different  preselected  directions  in  the  inertial  space.  The 
object  is  to  derive  the  equations  of  motion  capable  of 
describing  this  task. 

To  describe  the  motion  of  the  platform,  we  introduce  a  set 
of  inertial  axes  XYZ  and  a  set  of  body  axes  xyz  attached  to  the 
rigid  platform.  Then,  the  motion  of  the  platform  can  be 
defined  in  terms  of  three  translations  and  three  rotations  of  the 
body  axes  xyz  relative  to  the  inertial  axes  XYZ.  To  describe 
the  motion  of  the  flexible  appendages,  we  consider  a  typical 
appendage  hinged  at  point  e  and  regard  e  as  the  origin  of  a  set 
of  body  axes  xtyl2t  embedded  in  the  appendage  in  its 
undeformed  state.  Then,  the  motion  of  a  nominal  point  of  the 
appendage  consists  of  the  motion  of  xyz.  the  motion  of  xji^z, 
relative  to  xyz  and  the  elastic  motion  relative  to  xyy,  .  The 
system  and  the  various  reference  frames  are  shown  in  Fig.  1 . 

From  Fig.  1,  the  position  vector  of  a  point  in  the  rigid  body 
and  in  the  appendage  can  be  written  as 

8r  =  Bo  +  £  O'3) 

and 

8.  =  *=1-2 . -V  (l b) 

where  Ra  is  the  radius  vector  from  O  to  o,  r  is  the  position 
vector  of  a  nominal  point  in  the  rigid  body  relative  to  xyz, 
r„.  is  the  radius  vector  from  o  to  e.  r,  is  the  position  vector 
of  a  nominal  point  in  undeformed  appendage  relative  to 
xyy,  and  u,  is  the  elastic  displacement  of  that  point.  Vector 
R,  is  given  in  terms  of  components  along  XYZ  ,  r  and  r M 
in  terms  of  components  along  xyz,  and  r,  and  u ,  in  terms  of 
components  along  xyy,. 


and  that  of  a  point  in  the  typical  appendage  e  in  terms  of 
components  along  xyy,  is 

Y*  =  £.(&  +  x  £*)  +  (£*??  +  <ae)  *  (c,  +  u,)  +  ^ 

e  =  1 .2 . N  (46) 

where  aj,  is  the  angular  velocity  vector  of  axes  x  u  »  ,  E  is 
a  matrix  of  direction  cosines  between  the  xytzl  and  xyz  and 
v,  is  the  elastic  velocity  of  the  point  in  the  appendage  relative 
10  W.  -  &  =  •  In  the  maneuver  proposed,  the  angular 

velocity  vectors  <u ,  of  xyy,  relative  to  xyz  are  given,  so  that 
the  rotational  motions  of  the  appendages  relative  to  the 
platform  do  not  add  degrees  of  freedom.  The  only  degrees  of 
freedom  arise  from  the  rigid-body  translations  and  rotations 
of  the  platform  and  the  elastic  displacements  of  the 
appendages. 


The  equations  of  motion  can  be  obtained  by  means  of 
Lagrange's  equations  in  terms  of  quasi -coordinates  [Meirovitch 
(1988)). 


d  ,  dL  ,  -  dL 

dt  (  dVe  )  +  “  dVa 


(S  a) 


jL(1L)+  y 

dtK  do)  ’  0 


^  L  i  “Y  ^  L  f  r \T  *- 1  dL  ly  / » «i 
__  +  aj__(0  )  —  =Af  (56) 


<37",  „ 


e  =  1,2 . N  (5c) 


where 


L—T  —  V 


(6) 


The  velocity  vector  of  o  can  be  written  in  terms  of 
components  along  xyz  in  the  form 

Yo  =  CR0  (2a) 

where  C  is  the  matrix  of  direction  cosines  between  xyz  and 
XYZ  and  R,  is  the  velocity  vector  of  o  in  terms  of  components 
along  XYZ-  Matrix  C  depends  on  the  angular  displacements 
9,  (i  =»  1 ,2.3)  defining  the  orientation  of  axes  xyz  relative  to 
axes  XYZ.  Furthermore,  the  angular  velocity  vector  of  axes 
xyz  in  terms  of  components  along  xyz  is  given  by 


co  =  DO 


(26) 


where  9  is  a  vector  of  angular  velocities  9,  and  D  is  a  matrix 
depending  on  the  angular  displacements  0,  (i  =  1 ,2,3).  Figure 
2  shows  a  set  of  such  angular  displacements.  For  this  choice 
of  angles,  the  matrices  C  and  0  are  as  follows: 


C  = 


c9jc9  j 
-  cdySd) 
s9i 


cOjsO  j  +  rff|50jc0j 
c9\c9}  — 

-  s9\c9i 


sd\c9j  +  3 

Cd^d} 


(3a) 


s9  3 

0 

—  c9^s6j 

0 

s9i 

0 

1 

(36) 


where  c9,  =  cos  0j  and  s6,  =  sin  9,.  Note  that  this  choice  of 
angles  helps  us  avoid  singularities  at  the  initial  stage  of  the 
motion,  0,  =  0. 


In  view  of  the  above,  the  velocity  vector  of  a  point  in  the 
rigid  body  in  terms  of  components  along  xyz  is  simply 

Vr  =  Vg  +  co  x  r  (4a) 


is  the  Lagrangian  in  which  T  is  the  kinetic  energy  and  V  is  the 

A  A 

potential  energy,  L,  is  the  Lagrangian  density  and  T .is  the 
kinetic  energy  density,  both  for  appendage  e.  and  ££ ,  is  a 
matrix  of  homogenous  differential  stiffness  operators.  The 
terms  F  and  M  on  the  right  side  of  Eqs.  (5a)  and  (5b), 
respectively,  are  the  force  and  torque  vectors  on  the  platform, 

A 

both  in  terms  of  components  along  axes  xyz,  and  the  term  U. 
on  the  right  side  of  Eq.  (5c)  is  a  distributed  force  vector  on 
appendage  e  in  terms  of  components  along  xytzl.  Equations 
(5)  are  hybrid  in  the  sense  that  Eqs.  (5a)  and  (5b)  are  ordinary 
differential  equations  and  Eqs.  (5c)  are  partial  differential 
equations.  It  should  be  noted  that  the  tilde  over  a  symbol 
indicates  a  skew  symmetric  matrix  with  entries  corresponding 
to  the  components  of  the  associated  vector  [Meirovitch  and 
Kwak  (1988a)|.  For  the  system  of  Fig.  1,  we  write  the  kinetic 
energy 

f  N  r 

T  =  j  PrvJY4Dr  +  £  y  p.y!Y4D. 

Jor  ,-t  jd, 

N 

=  y  m<Yo  Yo  +  foil's  +  y  <st1ys  +  £ly 

+  yf  +  (Yo  +  'X“)T£.r(^«r9?«  +  [  p&dD') 

Jf>,  JD, 


+  mT  Ej  l, ajt  +  (E,q)  +  cy,)7  p,{r,  +  u,)v^D,l  (7) 

J0, 


where 


N 

”>t  =  mr  +  Yjm>  1 


mr 


PrdDr.  m. 
Or 


j  p4D,  (8<j,6,c) 

J0. 
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N 

s,  =  sr  +  +  eJS'E')  m 

«=*  1 

sr=  I  P/dDr.  S,=  )  pt(rf  +  U')dDt  (8 ej) 

JDr  J0, 

N 

4  =  4  +  +  E?!'F<  -  7i'F'S'F-'  -  F'lE>7o')  (**) 

«=i 

4  =  |  PS^dD,,  /,=  j  ne(re  4-  +  ue)TdDe  (8 ij) 

JD,  JD, 

in  which  p,  and  pf  are  mass  densities  and  D,  and  0,  are  the 
domain  of  the  rigid  platform  and  of  a  typical  appendage, 

respectively,  m,  is  the  total  mass  of  the  system.  S,  is  a  skew 
symmetric  matrix  of  first  moments  of  inertia  for  the  system  and 
/,  is  the  inertia  matrix.  For  simplicity,  we  assume  that  the 
potential  energy  is  due  entirely  to  elastic  elTects,  in  which  case 
il  ran  be  written  in  the  form 

■v  -v 

tZ  c  ^  ^ ]  =  ZJ  <9> 

M-lPs 


+  + 1 


N  . 

s,ya  +  its  +  Z  p,l'o.FJ  +  «.V.  +  ^)Vo,  =  s, 


-  ™!<°!  +  Z«  "  EI{2^'!>  ~  tr/«'J8)  +  2r0,£,#Ta®AI  +  ■?„) 

r=l 

+  2EI  |  Pei', f  +  ^ADe)FeV 

JD, 

-  Pet'c'E-Ju'  +  E^e(re  +  ue)]v/D, 

JD. 


+  roeE[(oj'S'Qj'  +  co,  +  Se„co,) 
~  £  {<*>,!,<&,  +  /,at,  +  /e„co,)}  +  M 

P,(Fe¥o  +  ( Eerle  +  (',  +  “e^efe  +  *,} 


=  Pei  -  \F,<S\  FeiF0  ~  'o,<H) 

+  2[v,  -  (r,  +  u,)co,  +  co ,(r,  +  u,))Eeq>  -  2cu.y, 

iV 

-  +  ",  +  "«)(E,  +  &)}  -  &eUe  +  {/. 


(He) 


where  [  ,  ]  denotes  an  energy  inner  product  [Meirovitch 
(1 980)|.  We  note  that,  in  deriving  Eq.  (9),  the  boundary 
conditions  were  fully  considered  (Meirovitch  (I980)|. 


III.  HYBRID  NONLINEAR  STATE  EQUATIONS  OF 
MOTION 

For  control  purposes,  or  for  mere  integration  of  the 
equations  of  motion,  it  is  convenient  to  work  with  state 
equations.  As  the  stale  vector,  we  consider 

J  =  ZbI  eT  «r  uj  ...  vj  a/  vj  v[  ...  sH  Y  (10) 

which  represents  a  unique  combination  of  inertial  coordinates, 
angular  displacements,  elastic  deformations,  translational 
velocities,  angular  velocities  and  elastic  velocities,  the  last  four 
being  in  terms  of  components  about  the  body  axes.  In  view  of 
the  definition  of  the  state  vector,  one  half  of  the  slate  equations 
consists  of  the  kinemalical  relations 

?0  =  cV0,  0  =  D''q>,  u,  =  V,  ,  e  =  1,2 . N  (lla,h,c) 

Inserting  Eqs.  (7)  and  (9)  into  Eqs.  (5)  and  considering  Eqs.  (2) 
and  (3),  we  obtain  the  hybrid  nonlinear  Lagrange's  equations 
for  the  other  half  of  the  state  equations 

n  e. 

+  sT'ii  +  yXr|  PSA^t  -  -  +  "Sr" 

«-i  io, 

n  p 

+  +  “  2"«  PSAD,  +  "«‘4"« 

.-I 


e  =  I.2.....V 


(11/) 


where  tr  denotes  a  trace  of  a  matrix  and 


Ee,  =  [  PS  AD,  (12a) 

4v  =  f  P&(\  +  +  (K  +  Wi VD.  (1 2b) 


IV.  THE  DISCRETIZED  NONLINEAR  STATE 
EQUATIONS  OF  MOTION 

The  equations  of  motion  are  hybrid,  in  the  sense  that  the 
equations  for  the  rigid-body  translations  and  rotations  of  the 
platform  are  ordinary  differential  equations  and  those  for  the 
elastic  motions  of  the  appendages  are  partial  differential 
equations.  Moreover,  because  of  the  maneuver  angular 
velocity  vector  a>,,  which  is  a  given  function  of  time,  they 
possess  time-dependent  coefficients.  Control  design  of  systems 
described  by  hybrid  equations  is  not  feasible,  so  that  we  wish 
to  discretize  the  partial  differential  equations  in  space,  leaving 
us  with  only  ordinary  differential  equations.  To  this  end,  we 
express  the  elastic  displacements  as  linear  combinations  of 


space-dependent  admissible  functions 
lime-dependent  generalized  coordinates,  or 

multiplied 

by 

“.(Cr  ')  =  ♦«(&>?•(')  •  «=  l*2- 

..,N 

(13) 

(Ud) 


where  <t>,  is  a  matrix  of  admissible  functions  and  q,  is  a 
vector  of  generalized  coordinates. 

The  Lagrangian  equations  in  terms  of  quasi-coordinates  for 
the  rigid  body  motions  of  the  platform  remain  in  the  form  of 
Eqs.  (5a)  and  (5b).  On  the  other  hand,  inserting  Eq.  (13)  into 
Eqs.  (5c),  we  obtain  the  ordinary  differential  equations  for  the 
discretized  elastic  motions 

-  ijr  -»•  ‘-u~"  "4> 

where 


0^04 
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■Vic ■" nj V AM t .*% CA  !  ,*69 

£.  =  ?.  05) 

<?«=[  p  —  1 ,2,...,N  (16) 

JD, 

are  corresponding  vectors  of  generalized  forces. 

The  Lagrangjan  remains  in  the  form  (6)  hut  the  kinetic 
energy  and  potential  energy  change,  indeed,  introducing  Eq. 
(13)  into  Eq.  (7).  we  obtain  the  discretized  kinetic  energy 

T  =  \  <",VT0V0  +  vis!*  +  y  «r/,« 

N  , 

+  ^l  Y +  y />J 

+  (&  +  rLwfzJtf'Q,  +  j  PfQ'dDrp') 

Jr> 


+  f  f.tv.iw, 

y  ~  r 

+  r,ro,e't>'  +  p[mt  pf[-t!?(]V^)P{ 

-I 


—  oS  +  a,lrQ  =  V  -*■  Tjr«E^(.Se<uj,  +  rjj,Seaje  +  Sfl,cu,) 

«-l 


4-  ai 


TEjl'<S,  +  (£,«  +  a»e)r  f  n,(r„  4-  [<! \qe~])<t>'dr),pe  I  (17) 

Many  of  the  quantities  in  Eq(l7)are  defined  hy  Eqs.  (S).  with 
the  exception  of 

~  f_  ~ 

•5,=  P,<%  +  Oe?J>*0e  (18n) 

J  n 


I,-\  Pt(r,  +  [^.<7,])(^  +  (IS*) 

Jp 

» 

Moreover,  inserting  Eq.  (13)  into  Fq.  (9),  the  discretized 
potential  energy  has  the  form 

N  N 

=  y£s«7’M»  O’) 

«=l  e-l 


-  Eelre<2*  + 

+  % 


r  ~ 

+  rVt>J[<D,q,]7rfD4Er)co 

JD. 

-  _  f  __  •  - 

+  4-  ^,<J>f(['D,7ja)(,  -  ai,r,  -  [‘t.p,] 
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(21*) 


/-V2 


rsj 

-  wtl<t>'qt~!i)dD'\r.,w  +  M,p,  +  2H,(m,)p, 

+  IK  +  //,(«,)  +  «.(£,«)  +  H.(Ue)k. 

~Q,~  +  I 

Jn 


JD. 


where 


where 


K  =  C  *>.■*' J 


(20) 


Following  the  same  procedure  as  used  earlier,  the  discretized 
nonlinear  state  equations  can  be  written  as  follows: 

N  N 

m,Xo  +  Sj<2  +  2^T  E,r([7,a),]  -  S„)E,<u  +  ^AtP* 

>-l  «-i 

N 

+  + 

nrtcu¥o  ~  a,St9} 

r— I 

s 

-  f  +  X  +  ®  Ae«)  (21a) 


r>~> 


S,t  +  +  £{£,r(2£,/,  -  tr/,5.)  +  2rME4r([5ecuf]  -  5,„) 


/-V/  r-W 


-  2£.r  f  P,(K  +  C0,<?J)C<t>^,]dDf)  E, 

h. 


+  f  pXCKqj]2r^Ot  +  <P[[E,wjE,r,q>  (21c) 

JP. 

•we  =  f  ptfQjn,  (22 a) 

♦,  =  f  (22*) 

5.  =  f  (22c) 

J0. 

/7,(a)=  f  ptfSt'df),  (22 </) 

JD. 

H'(a)=  f  ptffo'dD.  (22e) 

J/>. 

The  state  vector  is  redefined  as 

x  =  CR0r  0r  ?r  52r  -  ?i»  JKr  9/  £ir  Pj  -  £iv]r  (23) 

In  addition,  Eq.  (I  Ic)  is  replaced  by 

Pt~Pt  •  '=1.2 . (24) 
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Equations  (lla.b).  (24)  and  (21)  represent  the  discretized 
nonlinear  state  equations. 


V.  NUMERICAL  EXAMPLE 

The  effect  of  nonlinearity  on  the  system  response  is 
illustrated  by  means  of  a  spacecraft  consisting  of  a  rigid 
platform  with  a  single  membrane-type  flexible  appendage  (Fig. 
3).  The  maneuver  of  the  appendage  relative  to  the  platform 
was  carried  out  by  means  of  a  smoothed  bang-bang 
[Meirovitch  and  Quinn  (1987)|  for  the  angular  acceleration, 
where  the  smoothing  of  the  bang-bang  was  done  to  reduce  the 
clastic  deformations  of  the  appendage.  The  elastic  vibration 
of  the  appendage,  treated  as  a  circular  membrane  clamped  at 
r  =  a,  was  represented  by  ten  degrees  of  freedom  in  the 
z-direction,  i.e.,  by  ten  admissible  functions  in  the 
discretization(-in-space)  process,  so  that  the  matrix  <t>.  in  Eq. 
(13)  is  3  x  10,  or 


Two  cases  with  the  same  numerical  data  for  the  elastic 
appendage  but  with  different  inertia  terms  for  the  rigid  body 
are  tested.  The  data  for  the  elastic  appendage  is  as  follows: 

m,  =  0.303  slugs  .  S,  =  (0  0  1.4l5)r  slugs -ft 


/.= 


7.712  0.0  0.0 
0.0  7.712  0.0 
0.0  0.0  1.294 

.2 


slugs  •  ft 


p;  =  0.0l  slugs  I ft  ,  a- 1ft,  c  =  20  ft  Is 
^  =  (00  1.0 )T  ft 
For  the  rigid  body,  we  consider: 

Case  1 

mr  =  134.15  slugs  .  5,  =  0  slugs  •  ft 
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r  186.021 

0.0 

0.0  ' 

0 

0 

0 

0.0 

186.021 

0.0 

0 

<t>\ 

0 

*2  - 

0 

•p\0 

and  Case  2 

0.0 

0.0 

3S7.733 

slugs  •ft1 
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A)(ffoir) 


S>2  =  *1 


■Wit')  COS  9  <t> 4 
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4>i  J2(P n<x)  sin# 


<t>!)  2 

where 


■W2I')  cos  2 9 
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J'PPqT0) 
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*2  JiP22?)  sin29 


m,  =  21 .464  slugs  .  Sr  =  0  r/ugr  •  ft 


8.943  0.0  0.0 

0.0  8.943  0.0 

0.0  0.0  14.309 


slugs  • ft 1 


\fitPc  a 


*2  = 


As  seen  above,  the  rigid-body  model  of  Case  1  has  large 
inertias  relative  to  those  of  the  flexible  body.  On  the  other 
hand,  the  mass  moment  of  inertia  of  the  rigid  body  of  Case  2 
is  almost  the  same  as  that  of  the  flexible  body,  although  the 
mass  of  the  rigid  body  is  sufTicienliy  large.  This  is  due  to  the 
fact  that  the  mass  moment  of  inertia  of  the  rigid  body  is  about 
the  center  of  mass  and  the  mass  moment  of  inertia  of  the 
flexible  body  is  about  the  hinge,  which  is  far  removed  from  the 
mass  center  of  the  flexible  body. 


In  the  above,  and  <f1  are  recognized  as  axisymmelric  modes 
and  the  other  as  antisymmetric  modes,  respectively  [Meirovitch 
(1967)].  The  vector  q ,  is  ten-dimensional.  The  arguments  of 
the  Bessel  functions  “of  the  first  kind  can  be  obtained  from 
Pma  =  2.405,  /?0}a  =  5.520,  0„a  =  3.832.  Pl2a  =  7.016, 
/?2la  =  5.136  and  P^a  -  8.417.  The  mass  matrix,  Eq.  (22a), 
and  the  stiffness  matrix,  Eq.  (20),  arc  10  x  10  and  have  the 
block  -diagonal  form 

Me  =  I ,  K,  =  A 

where  /  is  tire  10  x  10  identity  matrix  and  A  is  a  diagonal 
matrix  with  the  diagonal  entries 

A(l.l)  =  -4</*Ola)2  •  A(2.2)  = 

a  a 


The  two  cases  are  compared  with  results  obtained  by 
Meirovitch  and  Kwak  (1988a)  using  linearized  state  equations. 
Figures  4  through  7  show  time  histories  of  the  rigid-body 
translations  and  rotation  and  the  elastic  displacements  of  the 
membrane  at  the  center  and  the  point  defined  by 
x,  =  0  and  r  =  1 .5  ft  for  Case  1 .  The  figures  contain  responses 
for  the  uncontrolled  nonlinear  system  and  linearized  system. 
Figures  8  through  1 1  show  the  responses  for  Case  2. 

The  elastic  displacements  at  the  center  and  at 
r  =  0.5  a,  9  =  90°  are  calculated  by  using  the  following 
formulas: 


(  1.08684  <7,  -  1.65807(72) 


A(3,3)  =  A(7,7)  = 


'r=0.5<7, 0-90"  —  7 —  (0.73006<7|  +0.29919^2 

v JPt 


A(4,4)  =  A(8,8)  = 

a 

A(5.5)  =  A(9,9)  =  4<02i'3>2 
a 

A(6,6)  =  A(10,I0)  =  4  (P^)2 
a 

where  c  =  sjTJp,  ,  in  which  Tm  is  the  membrane  tension  and 
p,  represents  the  mass  per  unit  area  of  the  membrane. 
Moreover,  the  other  matrices  given  by  Eqs.  (18)  and  (22)  are 
given  in  the  Appendix. 


-  1.06926<75  +  0.90574<76  +  1.15061  q7  -  0.35635<7|) 

A  discretized  nonlinear  state  equation  is  solved  by  using  1MSL 
routine  DIVPAG. 


VI.  SUMMARY  AND  CONCLUSION 

Slate  equations  of  motion  for  a  spacecraft  consisting  of  a 
rigid  body  and  a  given  number  of  flexible  appendages  are 
derived.  To  this  end,  Lagrange  s  equations  of  motion  in  terms 
of  quasi-coordinates  proved  to  be  most  effective.  In  general, 
the  resulting  hybrid  equations  are  nonlinear  and  time-varying. 


SIM  MECHANICS  PAN-AMEWCA  1969 


Appl  M*ch  ftov  19M  Supptomant 


In  addition,  the  equations  contain  persistent  disturbances  due 
to  inertial  loading.  Because  numerical  simulation  by  means 
of  the  hybrid  state  equations  is  not  feasible,  discretization  and 
truncation  are  carried  out.  The  discretized  state  equation  can 
be  used  for  the  maneuver  and  control  of  the  spacecraft,  such 
as  when  the  object  is  to  change  the  reorientation  of  the 
appendages  while  suppressing  their  vibration.  A  numerical 
example  of  a  spacecraft  with  a  membrane-like  antenna 
illustrates  the  differences  between  the  responses  obtained  by 
using  the  nonlinear  and  linear  state  equations. 
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APPENDIX 

For  the  given  configuration,  the  position  vector  of  a  nominal  point  in  the 
membrane  appendage  can  be  expressed  as 


r,  =  [  r  sin  5  r  cos  8  h  ] T  {^4  —  I ) 


Inserting  the  admissible  functions  given  by  Eqs.  (26).  together  with  Eq.  (A-l). 
into  Eqs.  (22),  we  obtain 


and 

W,<2)  -  -(of  *  °\)t 


where  /  is  10  x  10  identity  matnx.  Moreover, 

«.(J>  -  [0] 

f  P,  t  *,  ?,  i  JD,  =  [  0  ] 

n. 

where  we  note  that  the  above  null  matrices  are  10  *  10  and  3xJ0,  respectively. 
In  addition. 
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z  Angular  Displacement*  and  Velocities  of  the  Rigid  Platform  FIG.  3.  The  Spacecraft  with  a  Membrane  Ante. 
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FIG.  4.  Time  History  of  the  Rigid-Body  Translations  (Case  I) 


FIG.  6.  Elastic  Displacements  at  the  Center  (Case  1) 


FIG.  5.  Time  History  of  the  Rigid-Body  Rotations  (Case  1) 


FIG.  7.  Elastic  Displacements  at  the  r  =  0.5a  and  9  -  90*  (Case  1) 
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Abstract 

This  paper  is  concerned  with  the  proolem  of 
reorienting  the  line  of  sight  of  a  given  number  of 
flexible  antennas  In  a  spacecraft.  The  maneuver 
of  the  antennas  is  carried  out  according  to  a  min¬ 
imum-time  policy,  which  implies  bang-bang  control. 
Regarding  the  maneuver  angular  motion  of  the  an¬ 
tennas  as  known,  the  equations  of  motion  contain 
time-dependent  terms  In  the  form  of  coefficients 
and  persistent  disturbances.  The  control  of  the 
elastic  vibration  and  of  the  rigid-body  motions  of 
the  spacecraft  caused  by  the  maneuver  is  imple¬ 
mented  by  means  of  a  proportional -plus-integral 
control.  The  approach  Is  demonstrated  by  means  of 
a  numerical  example  in  whicn  a  spacecraft  consist¬ 
ing  of  a  rigid  platform  and  two  maneuveri  ng  ‘lexi- 
Ole  antennas  is  controlled. 

1.  Introduction 

Certain  space  missions  involve  the  reorien¬ 
tation  of  tne  line  of  signt  in  a  ‘lexto'e  soace- 
craft.  In  many  cases,  tne  line  of  signt  can  oe 
regarded  as  being  fixed  relative  to  tne  undeformed 
structure,  in  wnich  cases  the  reorientation  of  tne 
line  of  signt  implies  tne  reorientation  of  tne 
whole  spacecraft  (Ref,  1-5'.  Quite  often,  now- 
ever,  t re  lire  of  signt  coincides  wifi  a"  axis 
fixed  -r  a  small  component  of  tne  spacecraft,  sucn 
as  a  flexible  antenna.  In  sue n  cases,  it  may  oe 
more  advisaole  to  retarget  only  tne  flexible  an¬ 
tenna  and  not  the  entire  spacecraft.  Retargeting 
of  the  antennas  can  be  achieved  by  attaching  the 
antenna  to  tne  platform  trough  a  ni»-e  a5  to 
permit  o’/ot'"g  aoout  two  ortrjgo"al  axes.  Retar¬ 
geting  of  t"e  antennas,  instead  of  tne  wnole 
spacecraft,  is  a  virtual  necessity  when  there  are 
several  antennas  that  must  be  retargeted  simulta¬ 
neously,  and  the  line  of  sight  of  each  antenna 
must  be  reoriented  in  a  different  direction. 

Pefororre  5  considers  a  spacecraft  consisting 
of  a  r-gid  platform  with  a  number  of  flexible 
antennas  ;cig.  1',  and  tne  mission  is  to  retarget 
independently  tne  line  of  sight  of  each  antenna 
relative  to  the  inertial  space.  The  maneuvering 
strategy  developed  in  Ref.  6  consists  of  stabiliz- 
i "g  fe  o’ic'crm  'e’ative  to  the  ’"a^val  soace 
ang  reo'’ eng i "g  c-e  line  of  signt  of  each  antenna 
relative  to  tne  platform.  For  given  target  inac¬ 
tions  of  tne  antennas,  the  maneuvers  can  be  de¬ 
signed  as  if  tne  antennas  were  rigid.  Of  course, 
in  actuality  tne  antennas  are  flexible,  so  that 
the  maneuvers  are  lixely  to  cause  elastic  vibra¬ 
tion  of  the  antemas,  *mcn  in  turn  will  induce 


perturbations  t"e  oiitform.  hence,  tne  control 
tasx  amounts  to  simultaneous  staoi 1 i zation  of  tne 

platform  relative  to  the  inertial  space  and  vibra¬ 
tion  suppression  in  the  retargeting  antennas. 

The  mathematical  formulation  consists  of  a 
hybrid  set  of  equations  of  motion.  In  the  sense 
that  there  are  six  ordinary  differential  equations 
for  the  rigid-body  translations  and  rotations  of 
the  platform  and  partial  differential  equations 
for  the  elastic  motion  of  each  antenna.  Tha  equa¬ 
tions  of  motion  are  not  only  hybrid,  but  th*  ma¬ 
neuvering  of  the  antennas  relative  to  the  platform 
introduces  time-dependent  coefficients  Into  the 
equations.  Moreover,  the  equations  contain  terms 
-“fleeting  persistent  disturbances  caused  by  iner¬ 
tial  forces.  If  the  mass  of  the  antennas  is  small 
“elative  to  the  mass  of  tne  platform,  then  tne 
equations  of  motion  can  be  regarded  as  linear. 
Because  control  of  systems  governed  by  sets  of 
hyo“id  differential  equations  cannot  be  readily 
designed,  even  wnen  tne  equations  are  linear,  it 
!s  necessary  to  discretize  tne  partial  differen¬ 
tial  equations  in  space,  which  can  be  carried  out 
by  the  classical  Rayleign-Rltz  method  or  the 
finite  element  method  (Ref,  7).  Reference  6  pre¬ 
sents  tne  mathematical  formulation,  as  well  as  an 
example  illustrating  the  maneuvering  of  a  single 
flexible  antenna. 

This  paper  extends  the  work  of  Ref.  6  in 
several  respects.  In  the  first  place,  it  treats 
the  problem  of  retargeting  several  antennas  simul¬ 
taneously,  and  not  just  of  a  single  antenna.  In 
addition,  it  adaresses  the  problem  of  persistent 
d ’ st uroances  by  attempting  to  mitigate  tneir 
ecfect  during  tne  maneuver.  To  cope  witn  xnown 
disturbances,  disturbance-minimization  control  is 
effected.  The  retargeting  is  carried  out  open- 
loop  using  a  bang-bang  control  law.  This  implies 
that  the  inertial  forces  arising  from  the  maneuver 
angular  accelerations  are  almost  constant,  except 
for  a  sign  change  at  one  half  of  the  maneuver 
period.  lf  the  maneuver  is  not  very  fast  compared 
to  tne  lowest  natural  frequency  of  the  antennas, 
then  the  control  gains  can  oe  determined  by  ignor¬ 
ing  the  time-dependent  terms  in  the  coefficient 
matrices.  This  permits  the  use  of  proport ional - 
o’js-irtegral  feedbacx  control  for  disturbance 
accommoaat’on  ,Aef.  5,  9  and  10).  Of  course,  in 
t"e  computer  simulation  of  the  naneuver  and  con¬ 
trol,  the  full  time-varying  system  is  considered. 

This  paper  contains  the  procedure  for  design¬ 
ing  the  feedbacx  control  in  the  presence  of  dis¬ 
turbances.  Tne  prpportional-plus-integral  control 
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procedure  described  above  Is  demonstrated  by  means 

of  a  numerical  example  involving  a  soacecraet 
consisting  of  i  *•  i' i  .'r.'in  i-a  ~«Q  *"e<,-'e 
oeams,  eacn  w'tn  png  -•nge-:  to  c-e  o',  at'orm 
and  tne  otner  f-ee  v'ig.  2;,  wnere  tne  oeams  are 
originally  parallel  to  tne  z-axis  of  tie  plat¬ 
form.  Tne  maneuver  consists  of  slewing  eacn  of 
tne  oeams  tnrougn  a  45®  angle,  one  aoout  tne  x- 
axis  ana  the  otner  about  tne  y-axis  of  tne  plat¬ 
form. 

2.  Equations  of  notion 

We  consider  the  motion  of  a  spacecraft  con* 
sitting  of  a  rigid  platfons  with  a  given  numoer  of 
flexible  appendages  hinged  to  the  platform  as 
shown  In  Fig.  1.  The  rotation  of  each  individual 
appendage  relative  to  the  platform  is  given,  so 
that  the  motion  of  the  system  consists  of  the 
rigid  body  motions  of  the  platform  and  tha  elastic 
motion  of  tha  appandagas.  To  dascribt  the  motion, 
we  introduce  a  sat  of  Inertial  axes  XYZ,  a  set  of 
body  axes  xyz  attached  to  the  platform.  Then,  tne 
motion  of  the  rigid  body  can  be  deftned  in  terms 
of  three  translations  ana  three  rotations  of  axes 
xyz  relative  to  tne  inertial  axes  xyZ.  To 
descrioe  tne  motion  of  tne  flexible  appendages,  we 
consider  a  typical  appendage  nmgea  at  point  e  ana 
regard  e  as  the  origin  of  a  set  of  body  axes 
*QSeze  embedded  in  tne  appendage  m  jndeformed 

state.  rhen,  tne  motion  of  a  nominal  comt  of  t'e 
appendage  consists  of  tne  motion  of  <yz,  cne 
motion  of  *ejeze  relative  to  xyz  and  tne  elastic 

motion  relative  to  *~yaza.  The  system  and  tne 

various  reference  frames  are  snown  in  Fig.  1. 

c',om  Fig.  1,  tne  position  vector  of  a  point 
in  tne  rigid  body  and  in  tne  appendage  can  ce 
written  as  Rr  *  iQ  *  r,  *#  *  *0  *  £oe  *  Z9  ' 

(e  »  1,2,... ,N),  where  Rq  is  the  radius  vector 
from  0  to  o,  r  it  tne  position  vector  of  a  point 
in  tne  »-igia  body  naiative  to  xyz,  'S  tne 
radius  vector  from  o  to  e,  r  is  tne  position 

vector  of  a  nominal  point  in  undeformed  appendage 
relative  to  *e/eze  *nd  u^  Is  the  elastic  dis¬ 
placement  of  that  point.  Vector  Rq  is  given  In 
terms  of  components  along  XYZ,  r  and  r  In  terms 
of  components  along  xyz,  and  r^  and  ue  in  terms  of 
components  along  ^eze.  The  velocity  vector  of  o 
can  be  written  in  terms  of  components  along  xyz  in 

the  form  v  *  eft  ,  where  C  is  the  matrix  of  direc- 
'0  -o 

fon  cps'^es  between  xyz  and  XYZ  ana  ’s  tne 

velocity  vector  of  o  in  terms  of  components  along 
XYZ.  Furthermore,  the  angular  velocity  vector  of 
axes  xyz  in  terms  of  components  along  xyz  is  given 

by  j  *  0$ ,  where  |  is  a  vector  of  angular  veloci¬ 
ties  9 ^  and  D  is  a  matrix  depending  on  tne  angular 


displacements  (1  »  1,2,3).  In  view  of  the 

above,  t-e  velpp-ty  vector  pf  a  pom:  ,n  t-e 
POPy  m  ten^s  pf  ;;-cjnents  along  -5 

P1/  j_r  *  l  *  w  «  n  jnd  fiat  pf  a  point  in  tne 

typical  appendage  e  m  tenms  of  components  along 

Wo'2*  is  V  *  £  v  -  j  «  r  j  -  - 

e  -e  e  -a  ~  ~:eJ  "gd  de 

*  'Te  *  ~e'  *  ~e  e  *  1 12... .  ,N) ,  where  j  is  tne 
angj'ar  velocity  pf  axes  £a  -s  a  matnn  of 

direction  cosines  between  axes  *g/e2e  ana  xyz 

and  v  is  the  elastic  velocity  of  the  point  In  the 
appendage  relative  to  x^z,,  v#  -  &  .  in  the 
maneuver  proposed,  the  angular  velocity  vectors  u»# 
of  tgygZg  relative  to  xyz  are  given,  so  that  the 

rotational  motions  of  the  appendages  relative  to 
the  platfone  do  not  add  degrees  of  frtedoa. 

The  equations  of  motion  are  hybrid.  In  the 
sense  that  the  equations  for  the  rigid-body  trans¬ 
lations  and  rotations  of  the  platform  are  ordinary 
differential  equations  and  those  for  the  elastic 
motions  of  the  appendages  are  partial  differential 
equations.  Moreover,  because  of  tne  maneuver 
angular  velocities  they  possess  time-dependent 

coefficients.  Control  design  of  systems  described 
by  hypna  equations  is  not  feasible,  so  tnat  we 
w'sn  to  discretize  tne  system  in  space.  To  tms 
ena,  we  express  tne  elastic  displacements  as 
linear  combinations  of  space-dependent  admissible 
functions  multiplied  by  t ime-dependent  generalized 
coordinates,  or 

u  'r  ,t)  *  4  ( r  ]  q  ( t )  ,  e  *  1,2 . N  (1) 

-  ■'•c  C?  's,<? 

wnere  »  is  a  matrix  of  admissible  functions  and 
e 

je  is  a  vector  pf  generalized  coordinates. 

The  equations  of  motion  can  be  obtained  by 
means  of  Lagrange's  equations  in  terms  of  quasi- 
coonmates  ';ef.  11!.  X  linearized  version  of 
Sucn  equations  -as  den,vec|  gef<  g  for  tie  same 
matnematical  model  as  tne  one  considered  here. 

From  Ref.  5,  we  obtain  the  state  equations  of 
motion 


x(t) 
ere 


Mt)x(t)  +  BU)f(t)  *  0 ( t ) d ( t ) 


■re  x(t)  *  [RQ  3  3l  32*“3n  -o  ~  -l  52***-P^ 
a  state  vector,  in  wmch  s  is  a  symbolic  vector 


angular  displacements  of  the  platform, 
x  ge(e  »  1,2 . N)  and 


M*-l  = 


B(t] 
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0 
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M‘l(t)B*(ti 
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N*l(t) 

(3a) 
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are  coefficient  matrices,  in  which 
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is  a  mass  matrix. 
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Is  a  matrix  in  which 
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Moreover, 


1  0*  0”  ...  0N 


0  c  *  3  ...  3 

B*(t)  *  00  c2  ...  0 
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relates  the  discrete  force  vectors  to  the  modal  force  vectors,  in  wnicn 
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In  addition. 


c  =  -  *  e '  ~  e  1 '  fe  '~e2  1  -  V-Cen  i 
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Is  a  vector  of  disturbances  and 


Vi>  *  /„  oe#Il2#8dDe  ( 13n) 

''  J 

X'.ai  *  r^a  *  5^fj'»edDe  13o) 

Je 

In  wnic.n  a  ’s  a  vector  representing  j  or  j  ,  3 

-  6  ~  1?  r 

and  3e  are  mass  densities  of  eacn  oody,  0r  and  De 

are  t.ne  domain  of  tne  rigid  platform  and  of  a 
typical  appendage,  respectively,  Is  tne  total 

mass  of  the  system,  St  Is  a  skem  matrix  of  first 
moment  of  inertia  for  the  system  and  It  Is  the 

inertia  matrix.  For  simplicity,  we  assumed  that 
the  potential  energy  is  due  entirely  to  elastic 
effects.  In  which  case  it  could  be  written  In  the 
form 

11,1 
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plays  the  role  of  a  control  force  vector,  *nere  m 

#  * 

tne  latter  F_  and  M  are  actuator  *prce  anq 
-o  ~o 

torque  vectors  acting  on  tne  platform  »ng  ‘  ,  are 

actuator  forces  acting  on  aooendages  e  at  points  1. 
Other  quantities  entering  into  tne  aoove  matrices 
are  as  f ol lows : 


m  =  n  *  m  .  m  3  :  j 

t  r  e  =  ,  e  r  -3  e  r 


.  a  a ,  n 
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where  <e  is  the  stiffness  matrix  (Ref.  7). 

3.  Oisturbance-Accoweodatinq  Control 

The  maneuver  consists  of  retargeting  antennas 
so  as  to  point  in  given  directions  in  the  inertial 
soace.  3y  stabilizing  the  platform  in  an  inertial 
space,  the  task  reduces  to  reorienting  the 
antennas  relative  to  tne  platform.  For  a  minimum- 
time  maneuver,  the  control  law  is  bang-bang,  which 
implies  that  the  angular  acceleration  of  an 
antenna  relative  to  the  platform  is  constant,  with 
tne  sign  changing  at  half  the  maneuver  period. 
Ideally,  tne  maneuver  should  not  cause  elastic 
deformations  ’n  the  flexible  appendages.  This  is 
lisely  to  rec.-'e  a  long  maneuver  time,  which  is 
in  conflict  with  the  minimum-time  requirement.  To 
reduce  elastic  vibration,  a  smoothed  bang-bang  can 
be  used.  Still,  elastic  deformations  are  likely 
to  occur,  which  in  turn  implies  perturbation  of 
Me  olat'orm  -'em  3  fi <eo  position  in  the  inertial 
space. 
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The  motion  of  the  system  is  governed  by  Eq. 
(2).  The  system,  is  characterized  by  two  factors 
that  distinguish  it  from  most  commonly  encountered 
systems:  it  is  time-varying  and  it  is  subjected 
to  persistent  disturbances.  Both  factors  arise 
erom  the  retargeting  maneuver  angular  velocities 

w  .  angular  accelerations  j  and  the  matrices  E- 
—  6  —6  " 

of  direction  cosines  (e  *  1,2 . N),  all  quanti¬ 

ties  being  known  functions  of  time.  Clearly,  tne 
disturbance  term  0(tjd(t)  in  Eq.  (2)  depends  on 

Me  maneuver  oo' My.  In  tne  case  of  minimum-time 
maneuver,  tne  polity  's  oang-oang,  wmch  implies 
Mat  m8  maneuver  angular  acceleration  is  constant 
over  ootn  halves  of  tne  maneuver  period.  If  the 
maneuver  is  relatively  slow,  then  the  disturbance 
is  constant  over  both  nalves  of  the  maneuver 
oeriod.  In  tnis  case,  we  can  use  proportional - 
plus-integral  ‘eedoack  control. 
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Introaucing  the  notation 
3(tlu(t)  *  3 ( t ) f ( t )  -  3 ( t ) d ! t ) 
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Ea.  (2)  can  De  rewritten  as 


:  is) 
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(16) 


(254,6) 


x(t)  ■  A(t)x(t)  ♦  B(t)u(t) 

Assuming  that  O'v  s-d  \  vary  s’cwly,  ^ 
tnat  0(t'd(t)  ' s  almost  constant  Cum  "g  t~e 

control  interval,  we  can  write 


Gp  -  ,  fi1  •  6j  -  G^A 

Equation  ' 2 S '  '•eorosent^  -,ne  ootimal  control  law 
£or  tne  t ime- i nvar- ant  S/Stem  subjected  to  jnxnown 
constant  disturbances,  anO  is  snown  as  jropor- 
fonal-pl os-integral  control. 


u(t)  *  lit)  *  fd(t)  '  l7 

Introducing  a  new  state  vector  defined  oy 
2  *  C*1"  u1”]^.  £qs.  (16)  and  (17)  can  oe  combined 

Into  the  expended  state  equation 

j(t)  -  A(t)z(t)  ♦  B(t)fd(t) 

where 

A(t) 


A(t)  ! 

B(t) 

.  B(t)  - 

O 

o  • 

0 

I 

t  matrices. 

Note  that  If 

A( 

(18) 


(19a, 6) 


B(t)  are  a  controllable  pair,  then  A ( t )  and  3(t) 
are  also  a  controllable  pair  (Ref.  3). 


We  consider  an  optimal  control  polity  -n  tie 
e  that  f  (t)  minimizes  tne  performance  neasure 

J  *  i  zT(tf)Hz(tf)  *  j  f  C z”(t;Q!tlz't' 


The  optimal  control  law  is 
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where  G(t)  represents  a  control  gain  matrix  in 
which  K(t)  satisfies  the  matrix  Rlcattl  equation 

k  *  -  <A  -  AT<  -  3  *  <3R*13"<  ’2 

jsi ng  Eqs.  (16),  (17)  and  (21),  it  is  easy  to 
verify  that 


[n  general,  the  above  control  law  cannot  oe 
used  for  tne  type  of  proolem  considered  here, 
when  tne  maneuver  is  relatively  slow,  however,  so 
tnat  tne  matrices  A  and  B  are  nearly  constant!  tne 
control  law  (25)  can  oe  used  with  satisfactory 
results.  In  this  case,  the  control  must  be 
regarded  suboptlnel ,  but  close  to  being  optimal. 

We  recall  that,  to  reduce  vibration,  it  Is 
advisable  to  use  a  smoothed  bang-bang  for  the 
rigid-body  maneuver  of  the  appendages,  instead  of 
an  ideal  bang-bang. 

4.  Ibmrlcal  Exeaele 

The  mathematical  model  consists  of  a  rigid 
platform  and  two  flexible  beam,  each  one  with  one 
end  hinged  to  the  platform  and  the  other  end  free 
(Fig.  2),  where  the  beams  are  originally  parallel 
to  the  z-axis  of  tne  platform.  The  maneuver 
consists  of  slewing  eacn  of  the  beams  through  a 
45*  angle,  one  about  the  x-axis  and  the  other 
about  tne  y-axis  of  the  platform.  The  beams  are 
discretized  in  space  by  using  three  admissible 
‘unctions  for  eacn  component  of  displacement.  Six 
actuators  are  used  for  tne  rigid  platform  and 
three  actuators  for  each  displacement  component  of 
both  beams.  The  latter  actuators  are  located  at  4 
ft,  7  ft  and  10  ft  from  tne  pivot  point,  the  third 
coinciding  with  the  tip  of  tne  beam.  Figure  3 
shows  the  maneuver  time  histories,  in  which  the 
angular  acceleration  represents  a  modified  bang- 
bang.  Figures  4a  a-d  4p  display  both  the  un¬ 
controlled  and  cont-olled  translational  and 
angular  displacements  of  the  platform,  respect¬ 
ively,  and  Figs.  5a  and  5b  show  the  tip  dis¬ 
placement  of  the  two  beams.  As  can  be  verified, 
the  maneuver  and  control  of  the  spacecraft  is 
suite  satisfactory.  *he  disturbance-accommodating 
cont-ol  is  earned  out  by  tne  proportional-plus- 
integral  control  approach. 


df  t  t<  In  obtaining  the  numerical  results,  the 

^  ■  (Gj  -  GjB  A)x  ♦  G,B  x  (23)  following  data  was  used: 


where  Gj  and  G-j  represents  submatrices  of  G 
correspond  ng  to  <  and  u,  correspondingly ,  and 

t  r  I  T 

8  *  (3  8 ) ~ *■  8 1  is  the  pseudo-inverse  of  B. 

Integrating  Eq.  (23),  we  obtain  the  optimal 
control  law 

f't:  *  ‘'0)  r  /  (G,  -  3,3fA)x  ot 

0  1 

*  '*  G,Bf*  dt  (24) 

'  0  " 

If  tf  is  Sufficiently  large  and  A  and  8  can  be 

assumed  to  be  constant,  then  the  gain  matrix  is 
constant.  Moreover,  if  x(0)  and  f(0)  are  zero, 
then  Eq.  (24)  yields 

t 

f ( t )  *  G  x  ♦  G .  f  x  dt  (25 ) 

p~  i  •  Q  - 


mr  *  134.15  slugs  ,  mg  *  0.1373  slugs 
Sr  =*  [0  0  or  slugs. ft  , 

Se  *  [0  0  0.9365 JT  slugs,  ft 

""196.021  0  0 

0  126.021  0  slugs,  ft2 

0  0  357 . 7  3  3_ 

6.243  0  0 

0  8.243  0  slugs. ft2 

0  0  0 


!e  5 


1  .  12  »  10  ft,  £1  »  3023.9  lb. ft 

rQl  *  [0  -1.0  0.5jT  ft 


where 


£q2  -  [0  1.0  0.5]T  ft 
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«nere  !  ts  tie  Identity  matrix.  Qf  course, 
consistent  witn  a  steady-state  solution  of  tie 

matrix  Riccatl  equation,  h  was  taken  as  zero. 
Finally,  for  control  design  purposes,  the 
coefficient  matrices  A  and  3  were  taken  as 
constant  and  corresponding  to  the  premaneuver 
configuration  of  the  spacecraft.  Of  course.  In 
implementing  the  control,  the  time-varying 
matrices  A(t)  and  8(t)  were  used. 
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5.  Sa— nr  and  Conclusions 

Certain  space  missions  involve  the  reorient¬ 
ation  In  an  Inertial  space  of  the  line  of  sight  of 
certain  small  flexible  components  of  a  spacecraft, 
such  as  flexible  antennas.  In  sucn  cases,  '■egara- 
ing  the  main  spacecraft  as  a  rig-.-  platform,  a 
sensible  strategy  is  to  stabilize  the  platform 
relative  to  trie  inertial  space  and  '■eorient  tie 
line  of  sight  of  the  various  antennas  relative  to 
the  platform. 

’’his  paper  is  concerned  with  toe  problem  of 
retargeting  several  antennas  simultaneously,  wnile 
suppressing  any  rigid-body  and  elastic  perturbat¬ 
ions  caused  by  the  retargeting  maneuver.  rhe 
retargeting  was  carried  out  open-loop  using  a 
smoothed  bang-bang  control  law.  'his  implies  that 
tne  inertial  forces  arising  f~on  the  maneuver 
angular  accel erati ons  are  a’most  constant,  except 
for  a  sign  change  at  one  naif  of  t.ne  maneuver 
period.  If  the  maneuver  is  not  very  fast  compared 
to  the  lowest  natural  frequency  of  the  nonmaneu¬ 
vering  antennas,  then  the  control  gains  can  be 
determined  ay  Ignoring  tne  t i ie-deoende"t  terms  in 
;ne  coefficient  matrices,  '-'is  cermitted  the  use 
of  oroport iona l -p 1  us- 1 ntegra 1  ‘eedbacx  control  ‘or 
disturbance  accommodation.  Of  course,  in  the 
computer  simulation  of  the  maneuver  and  control , 
the  full  time-varying  system  was  considered. 

A  numerical  example,  in  whicn  a  spacecraft 
consisting  of  a  rigid  olatform  and  two  flexible 
antennas  jngongoing  naorientation  in  different 
planes  was  controlled,  demonstrates  the  approach. 
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